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■ Abstract. We investigate the impact of three-body forces on the transverse momentum distribution of 

partons inside the proton. This is achieved by considering the three body problem in a class of hyper- 
central quark potential models. Solving the corresponding Schrodinger equation we determine the quark 
wave function in the proton and with appropriate transformations and projections we find the transverse 
momentum distribution of a single quark. In each case the parameters of the quark potentials are adjusted in 
order to sufficiently describe observable properties of the proton. Using a factorization ansatz we incorporate 
the obtained transverse momentum distribution in a perturbative QCD scheme for the calculation of the 
cross section for prompt photon production in pp collisions. A large set of experimental data is fitted 
^ C"| using as a single free parameter the mean partonic transverse momentum. The dependence of (kr) on the 

O |' collision characteristics (initial energy and transverse momentum of the final photon) is much smoother 

when compared with similar results found in the literature using a Gaussian distribution for the partonic 
transverse momenta. Within the considered class of hyper-central quark potentials the one with the weaker 
dependence on the hyper-radius is preferred for the description of the data since it leads to the smoothest 
mean partonic transverse momentum profile. We have repeated all the calculations using a two-body 
potential of the same form as the optimal (within the considered class) hyper-central potential in order to 
check if the presence of three body forces is supported by the experimental data. Our analysis indicates 
that three-body forces influence significantly the form of the parton transverse momentum distribution 
and consequently lead to an improved description of the considered data. 

ON 

1 PACS. 12.39.Jh , 12.39.Pn, 12.38.Qk, 12.39.Ba, 13.85.Qk 

^6 _ 

|Q | 1 Introduction relativistic inter-quark potential 

V(r) = A + Br 01 (1) 

The description of the hadron spectra is still an open ques- 
tion in theoretical physics. To date, the most important w i t h A, B appropriate constants, has been successfully 
^ . progress in this direction is based on lattice QCD, QCD usec } m the literature for the description of the heavy 
H i sum rules and potential models. Despite being more fun- quar k me son wave function as well as the clearly relativis- 
damental, lattice QCD and QCD sum rules can only lead tic ss states 1 . The same model was used later in [2] in 
to rough description of hadronic spectra. Potential mod- or Her to obtain baryonic spectra with very good results, 
els on the other hand, albeit not fundamental, have been with the progress of lattice QCD it became possible to 
proved to be very successful even in the non-relativistic determine effective inter-quark potentials based on hrst 
approximation. principles 0E]. D uring the late eighties it was realized 
Since the late seventies several attempts have been that genuine many body interactions could also play an 
made in this field using different forms for the quark-quark important role in the determination of baryon properties 
interacting potential, leading in many cases to a very good [7] . To this end it has been proved to be very efficient 
description of the hadronic states. The main ingredient in to express the interaction between the quarks in terms of 
all these models is the presence of a confining part in the hyper-radial potentials [51r0l fT0lfTT) . In this treatment V(£) 
inter-quark potential, which is independent of flavor and is in general a three body potential since the hyper-radius 
spin. £ depends on the coordinates of all three particles. Such 
Although in light baryons, as the proton, relativistic potentials have been extensively used for a consistent de- 
effects are expected to play a role, there are several non- scription of a large set of hadronic observables which be- 
relativistic treatments leading to satisfactory description sides their spectra include the photo-couplings (T3], the 
of the proton properties In particular the non- electromagnetic form factors and the strong decay ampli- 
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tudes [T3"lll4] . Recently it has been proposed that within 
the potential model approach one could also obtain an es- 
timation of the transverse momentum distribution of par- 
tons inside the hadrons g(kr) [IS]- The resulting prob- 
ability density, characterized by a non-Gaussian shape, 
was then used within the framework of perturbativc QCD 
for a phenomenological description of the cross section for 
7T production in pp collisions. Interestingly enough, this 
treatment turned out to be very efficient in the description 
of the experimental data resolving several unsatisfactory 
issues present in the usual approach involving a Gaussian 
form for gikx). In |15) the treatment was based on a two- 
body potential of the form ([T]) while in a later work the 
MIT bag model [TB] was used in a similar manner to ob- 
tain g(kx) and subsequently to describe successfully the 
prompt photon production in pp collisions |17] . The results 
of these two works indicate that confinement, asymptotic 
freedom and/or relativistic description make an imprint 
on the intrinsic transverse momentum distribution of the 
constituent quarks, detectable in the cross section of pp 
collisions. As mentioned in [17] this fact could give an 
explanation for the systematic discrepancy between theo- 
retical next-to-leading order (NLO) calculations [TS\ and 
experimental data [19,20] of inclusive single photon pro- 
duction (pp — > jX). The observed gap for this process 
is particularly significant in fixed target experiments and 
cannot be filled even after taking into consideration cer- 
tain large contributions to the partonic hard scattering 
cross section to all orders in perturbation theory, using 
the threshold resummation technique [2"Tll2"2"inZ3"] . 

In the present work our main interest is to explore 
if genuine three-body effects may also influence the trans- 
verse momentum distribution of the quarks inside the pro- 
ton in a way that it is detectable in cross section data. 
Therefore, we initially attempt a consistent description of 
the ground state wave function of the proton as a three 
quark bound state within a class of hyper-central quark 
potential models of the form: 



Our work is organized as follows. In section 2 we in- 
troduce the hyper-central description of the three body 
problem in the non-relativistic case using the potential 
V(£) — Ak + £>fc£ fc ■ We solve numerically the hyper-radial 
Schrodinger equation for the 4 different values of k men- 
tioned above. In each case we determine the parameters 
Ak and B k having as criterion the exact description of 
the ground and first excited state energy of the proton. In 
section 3 we calculate the intrinsic transverse momentum 
distribution g(k,T , (^t)) of partons inside the proton for 
different k. In section 4 we present the numerical results 
from the best fit of the pp — > cross section data at vari- 
ous energies and transverse momenta of the produced pho- 
ton, using the single parameter distribution g(kx, (&t)) 
within the framework of perturbative QCD. Finally, sec- 
tion 5 contains our concluding remarks. 



2 Hyper-central potential for the proton 

We start our study by considering the proton as a bound 
state of three constituent quarks. After fixing the center 
of mass, the three particle configuration is described by 
the Jacobi coordinates: 



€ 2 = 



£ 1 = x 2 - xt 

2X3 — X\ — X2 

7! 



(3) 



Instead of Xi (i = 1,2,3) one can introduce the hyper- 
spherical coordinates, which are given by the angles fl\ = 
(6i,(pi) and fli = (02,4>2) (Qi, <pi are polar and azimuthal 
angles of vector i = 1,2) along with the hyper-radius 
£ and the hyper-angle x> defined by the relations: 



6 



£ cosx 
: £ sin x 



(4) 



V(0 =A k + B k e (2) 

where A kl B k are constants, k G R + and £ is the hyper- 
radius. Due to the many body nature of the problem in 
the general case the wave functions can only be obtained 
numerically One interesting exception is when the quark- 
quark confining potential is harmonic, allowing for analyt- 
ical solutions. In addition, the harmonic model supplies a 
convenient classification scheme of the baryon resonances 
in terms of shells [3] . The parameters A k and B k are cho- 
sen in order to fit the proton's ground and first excited 
state energy. To ensure consistency we also estimate the 
proton's charge radius. Having fixed Ak, B k we determine 
g(kx) for several values of k (k = 0.1,0.5, 1,2) and then 
we use the standard treatment within perturbative QCD 
for the calculation of the cross section for prompt photon 
production in a pp collision experiment. We compare our 
results with those of [T51[T7] in order to extract informa- 
tion concerning the presence or not of traces of three-body 
effects in g(kr) traceable through the considered experi- 
mental data. 



In this model we consider three identical quarks of mass 
to. Then the Hamiltonian can be written as: 



H 



2/j. 



+ v(0 



(5) 



where fj, = % and the potential depends only on £ (hyper- 
central). Since £ 2 = £1 + £| the interaction in Eq. [S]is not 
a purely two body interaction but contains three body 
contributions. The presence of three quark forces could 
be suggested by the existence of a direct three gluon in- 
teraction which is one of the fundamental features of the 
non-abelian nature of QCD. In fact all these many body 
terms can be included, effectively, in an appropriate hyper- 
central potential [TT] . 

Using hyper-spherical coordinates the Hamiltonian of 
the three body problem can be written as: 



H = -±- 

2fi 



1 d 



e at v at 







v(0 

(6) 
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where: 



; 2 X d X 



■ 2 „ d 
sin 2x^- 



X 



sin 2 x 



(7) 

in which J is the angular momentum of the subsystem of 
particles 1, 2 and L is the angular momentum of particle 3 
with respect to the center of mass of the two body subsys- 
tem (1, 2). £ 2 (J?i, Q 2 , x) is t ne Casimir operator of the six 
dimensional rotation group 0(6) and its eigenfunctions are 
the hyper-spherical harmonics Y\j mj i mt (f2i, f2 2 ,x)- That 
is: 



C 2 (r2 1 ,f2 2 ,x)Y\j mj i mi (f2i,n 2 ,x) = 
= -\{\ + 4)Yx jmj i mi {Oi,n2,x) 



(8) 



where the grand-angular quantum number A is given by 
A = 2n + j + l, n = 0, 1, 2, . . . and j, I are the angular mo- 
menta associated with the J and L operators respectively. 
The solution of the Schrodinger equation in this case has 
the following form: 

& N\jmjlmi{£,X> di,(pi,0 2 , <t>2) = 

— N]y\j mj i mi R]si\(!;)Y\j mj i rni (Qi,(f)i,Q2,4>2,x) (9) 
where 

*A im . /mi (0i - <h,h,<h,x) = cos J+ ^ X sin^ xP l n +hj+i ■ 
■ (cos2 X )Yp(9 1 ,ct >1 )Y l ^(e 2 , ( p 2 ) (10) 

i+- i+- 

and V n 2 ' 2 (cos 2x) are the Jacobi polynomials. For the 
hyper-radial many-body interaction between equal mass 
particles we use the general form @ leading to the radial 
equation: 



H £ 5 <9£ V d£ 



Rn\(C = o 



(11) 



Introducing new variables 



P = V^T" > £o = ( — ) , e = -Bjva - A k 



lh W = 2i^i 1 



2/^e 2 

and the new function: 

un\(p) = p 5/2 Rn\{p) 
As a result Eq. (fTTj) becomes: 



(12) 



(13) 



rf2 (A+i)(A + i)_ pfe+Tye 



w WA (p)=0 (14) 



determining the energy eigenvalue problem to be solved. 
As discussed previously the non trivial part of the calcu- 
lation is the solution of the radial equation (Tl4^ . which in 
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Fig. 1. The ground state wave function of the proton using 
the hyper-radial potential ^ for four different choices of k: 
0.1 (solid line), 0.5 (slashed-dotted line), 1.0 (dotted line), 2.0 
(slashed line) for the potential V(p) = a k + f3 k p k '■ The thick 
solid line correspond to the MIT bag model. The inset displays 
the ground state wave function for k — 0.1 in an appropriate 
scale. 



the general case (arbitrary k) can be obtained only numer- 
ically using the Numerov algorithm. Actually, solving (14) 
we determine simultaneously the eigenvalues do = Wea 
and (Ti = Wei and using the equations: 



r(fi, B k )A k = r(p, B k )Ei - 



2ai 
E, 



i = 0, 1 ; T(tx,B k ) = 2^ 



(15) 



where E Q = 0.938 GeV and E 1 = 1.440 GeV are the pro- 
ton's ground and first excited state energies respectively, 
we find the constants A k and B k of the hyper-radial po- 
tential in ([2]). These quantities are necessary in order to 
get estimations of dimensionfull observables from our cal- 
culations. In Figs. II 121 we present the wave functions of 
the ground state and the first excited state respectively 
(properly normalized for dimensional reasons) for the four 
different values of k mentioned above. The insets display 
in more detail, using a suitable scale, the form of the 
ground state (Fig. [T]) and the first excited state (Fig. [2]) for 
k = 0.1. Using the ground state wave functions shown in 

_ 1 C 2 1 

Fig. [T] we can calculate the hyper-radius (£ 2 ) 5 = ^-(p 2 ) 5 

of the proton ((p 2 )^ being the corresponding dimension- 
less quantity). 

The hyper-radius can be related with the experimen- 
tally accessible charge radius r c h = \J (x 2 ) , i = 1,2,3 
(assuming (xf) = (a; 2 ) = (xj)) through: ^/{xf) ' ch = 

\ \J (£ 2 ) ■ I R Table I we summarize the results for the pro- 
ton charge radius using the four different values of k pre- 
viously mentioned. 

According to Table I, based on the value of the proton's 
charge radius no distinction between the four considered 
cases of k is possible. 
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Fig. 2. The first excited state wave function of the proton 
using the hyper-radial potential ((2)1 for four different choices of 
k: 0.1 (solid line), 0.5 (slashed-dotted line), 1.0 (dotted line), 
2.0 (slashed line) for the potential V(p) = a k +/3kp k ■ The thick 
solid line correspond to the MIT bag model. The inset displays 
the ground state wave function for k — 0.1 in an appropriate 
scale. 

Table 1. The proton charge radius for different choices of the 
exponent k in the hyper-radial potential. 



k 


r c h (fm) 


0.1 


0.60343 


0.5 


0.60695 


1 


0.60901 


2 


0.61002 



3 Intrinsic Transverse Momentum 
Distribution g(k T ,{kT)) 

As a next step we determine, for each choice of k, the 
single particle transverse momentum distribution in the 
ground state. In order to proceed we first have to calculate 
the ground state wavefunction in the momentum space 
(conjugate to the space £j , i = 1, 2, 3) as: 



0(fc c ) = N J d^dxdf2 1 df2 2 ^ 5 cos 2 x sin 2 \ • 

■ exp[-«fc 5 £ cos x cos 6x]W N \j mj i mi (£) (16) 

where &N\jm«lm l is the full eigenfunction. A convenient 
representation of the integrals in eq. (fTB|) is achieved in 
the reference frame where = (0, 0, k^, 0, 0, 0). 

It is useful to determine the transformation of the mo- 
menta fc^ to the Cartesian momenta kf. 



^2 



1 

1 



1 



(fei - k 2 



(fei + fe 2 - 2fc 3 ) 



(17) 



Fig. 3. Transverse momentum distribution of a parton inside 
the proton corresponding to the different values of k in the 
hyper-radial potential model @. The thick solid line corre- 
sponds to the MIT bag model. 



The above expressions is simplified, in the center of mass 
frame where k^ 3 — and: 



kl 



kf 



ki ■ k 2 



The two particle density p{k\, k 2 ) is then given by: 



p{k 1 ,k 2 ) = 4>(kt) 



kf + kl + k 1 ■ k 2 ] 



(18) 



(19) 



Finally from Eq. [19] we obtain the one particle transverse 
momentum density g(kx) as: 



r+oo 

g(k T ) = 4vr / dk, 



+ 00 

2 



(20) 



with z = cos uj where cj is the angle between fei and 
k 2 . The integration in equations [16] and [20] can be per- 
formed to a great accuracy using a mixture of Gauss- 
Kronrod quadrature and VEGAS Monte-Carlo integration 
algorithm. In Fig. [3] we present the intrinsic transverse 
momentum distribution l-K^gikr) of a parton inside the 
proton obtained from the ground state wave function cor- 
responding to each of the four different choices of the ex- 
ponent k in @ . As expected for increasing k the maximum 
of the distribution becomes broader while the tail tends 
to be more abrupt. 

One possibility for testing the phenomenological rele- 
vance of three-body forces is to consider their influence in 
the description of physical processes, like prompt photon 
production in pp-collisions, within the framework of per- 
turbative QCD. Our strategy is to find the value of k in @ 
for which we achieve the best description of experimental 
data and then to compare our results with those obtained 
using g(kr) determined through the MIT bag model [17] 
or the corresponding two-body potential |15j . 
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4 Numerical Results 

We consider prompt photon production in pp-collisions 
as the appropriate process for checking the influence of 
three-body partonic interactions in the phenomenology of 
proton collisions. In fact this process is optimal for this 
purpose since it is not affected by experimental ambigui- 
ties caused by final state hadronic interactions. In order to 
proceed we use the phenomenological scheme proposed in 
p~5llT7] for the calculation of the differential cross-section 
for inclusive 7-production. This scheme incorporates par- 
tonic subprocesses according to perturbative QCD, par- 
tonic effects in the proton described through the longi- 
tudinal parton distribution functions (PDF) and effects 
due to the intrinsic transverse momenta of the partons 
described through g(kr). A simplified phenomenological 
approach is adopted, in which it is assumed a factoriza- 
tion between longitudinal and transverse momentum par- 
ton distributions [2~4"I|25| . Although such an assumption 
seems reasonable from a statistical point of view, since 
the longitudinal momenta of the partons may differ by 
orders of magnitudes from the corresponding transverse 
ones, its validity based on first principles remains under 
question [33ll34l . Despite this fact this factorization asatz 
turned out to work sufficiently well in the case of cross 
section calculations for prompt photon production using 
the MIT bag model [T7]. Here, as we are interested in 
comparing results obtained using different inter-partonic 
interactions for the description of the proton wave func- 
tion, it is necessary to use exactly the same treatment as 
that introduced in [T51[T7] . The calculations are performed 
in next-to-leading order (NLO) of perturbative QCD and 
the cross section for single photon production is given by: 

E 1 -p-(pp ->• 7 + X) = K(p T , Vs) ^ J dx a dx b ■ 

abc 

■fa/p(x a , Q 2 )fb/ P (xb, Q 2 ) ■ -> cry)6(s + i + u) (21) 

7T dt 

where f i/p (i = a, b) are the MRST2006 NNLO longitudi- 
nal parton distribution functions (PDF) for the colliding 
partons a and b as a function of longitudinal momentum 
fraction Xi and factorization scale Q [26]. ^? is the cross 
section for the partonic subprocesses as a function of the 
Mandelstam variables s, t, u [35] ■ The higher order correc- 
tions in the partonic subprocesses are effectively included 
in (|2ip through the if-factor, appearing in the right hand 
side, which depends on the transverse momentum of the 
outcoming photon and the beam energy [27] . 

At this point we should mention that although part of 
the fey-effects is unavoidably included in the NLO calcu- 
lations, here we are studying the non-perturbative origin 
of such effects. To this purpose we are using a minimal 
modification to the standard approach in order to obtain 
an upper bound for such non-perturbative effects. Follow- 
ing this reasoning we attempt to describe experimental 
data introducing partonic transverse degrees of freedom 
through the replacement [231128) : 

dxi fi/ p (x l ,Q 2 ) — ► dx i d 2 k T . i g(k Tii )f i / p (x i ,Q 2 ) (22) 



in the PDF of the colliding partons (i — a,b). To avoid 
singularities in the partonic subprocesses we introduce a 
regularizing parton mass [291130] with value close to the 
constituent quark mass m q = 0.3 GeV in the Mandelstam 
variables appearing in the denominator of the correspond- 
ing matrix elements. In fact m can be chosen in the range 
[0.1, 1.0] GeV without affecting the following analysis. Us- 
ing the distribution g(kr) obtained in the last section it 
is straightforward to calculate the cross section (|2ip . 

We start our numerical investigations calculating the 
differential cross section for PHENIX data [21] on prompt 
photon production with transverse momentum pr at RHIC 
(y/s = 200 GeV). At this step of the analysis we use 
the PHENIX data since, due to the very high beam en- 
ergy, the ^7= ratio is expected to become very small in- 
dicating the presence of non-perturbative QCD processes 
where fey-effects are expected to be relevant. We perform 
four sets of runs, each one using a different distribution 
g(kr) (see Fig. [3]) associated with the different values of k 
(k = 0.1, 0.5, 1, 2) in the potential ©. Varying the mean 
transverse momentum (fey) we fit in each case all the avail- 
able PHENIX data for different pr of the produced pho- 
ton. As a result an one-to-one relation of (fey) with pr, 
for each <?(fcy) used, is established. In Fig. [4] we display 
graphically this relation for the four considered cases. In 
general the variations of the (fey) dependence on pr are 
not as large as when an ad-hoc Gaussian g(kr) is used 
(El)- However, in order to achieve a comparison between 
the different models we impose the following two criteria: 

— compatibility of the fitted (fey)-values with the geo- 
metrical properties of the proton, and 

— smoothness of the relation between (hr) and pr 

To make the first requirement more quantitative we cal- 
culate the pr averaged uncertainty <X(fc T ) = A(hr) of (hr) 
for each of the considered models and we compare it with 
tt^— obtained using the proton charge radius shown in ta- 
ble I. Assuming that <J(k T ) describes successfully the fey- 
fluctuations then the ratio: 

Rg = ^f- (23) 
2r ch 

can be used as a measure of the consistent description of 
the proton's geometry (size) within the considered model. 
Optimally we expect Rg ~ 1 while deviations may origi- 
nate from the type of the inter-quark potential, the pres- 
ence or not of three-body forces and the relevance or not of 
relativistic effects. The second criterion is quantified intro- 
ducing a non-smoothness parameter Rns defined as the 
average slope variation squared in adjacent py-intervals. 
To be more precise one uses a linear approximation for the 
function (fcy)(py), as determined by the pairs (pr, (hr)), 
found through the fitting of the experimentally observed 
cross section for each considered model, to estimate the 
slope Si in the i-th py-interval. Then Rns is given by: 

1 N-l 

Rns = j^Y,^~ s ^ 2 ( 24 ) 

i=2 
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From this definition it is clear that with increasing Rns 
the associated function (fcx)(pT) becomes less and less 
smooth. The results for the quantities Rg and Rns, cal- 
culated using g(kx) obtained from the four different po- 
tential models discussed above, are summarized in Table 
II. For comparison we also include in this table the values 
of Rg and Rns found using g(kx) determined by the MIT 
bag model [17] as well as by solving the three body prob- 
lem with two-body interactions of the form (|TJ) . In order 
to be complete we give the values of these quantities found 
in the case of using the usual Gaussian g{kx) in the cross 
section calculations. 



Table 2. 

els. 



The quantities Rg and Rns for the considered mod- 



Model 


Rg 


Rns 


k = 0.1 (3-body) 


1.27 


0.17 


k = 0.5 (3-body) 


1.45 


0.18 


k = 1 (3-body) 


1.67 


0.30 


k = 2 (3-body) 


1.88 


0.57 


k = 0.1 (2-body) 


2.97 


1.18 


MIT bag 


0.85 


0.04 


Gaussian 


4.06 


3.65 



According to Table II it is evident that the hyper-radial 
potential with k = 0.1 leads to more consistent values for 
Rg and Rns than the other three choices (k = 0.5, 1.0 
and 2.0). Clearly for g(kr) obtained from © with k = 0.1 
the fluctuations of (kx) are smaller and closer to the ex- 
pectations for the proton size based on Heisenberg's un- 
certainty relation (Akx ~ 0.33GeV). Comparing the re- 
sults for the k = 0.1 hyper-radial potential with those for 
the similar 2-body potential we conclude that three body 
forces are important for a consistent description of the 
partonic transverse momentum effects in the proton. In 
addition, the partonic transverse momentum distribution 
determined using the MIT bag model leads to the best val- 
ues for Rg and Rns suggesting that relativistic effects also 
influence significantly the transverse momentum structure 
of the proton. Finally, it is important to trace the behavior 
of each model to the characteristics of the corresponding 
distribution gikr). As it is clearly seen from Fig. [5] the 
best description, of the experimental data, is achieved us- 
ing the model resulting in the greatest variance of kx for 
given (kr) (MIT bag [17 ). It is also interesting to notice 
that this characteristic depends smoothly on the exponent 
k of the potential, becoming more pronounced for lower 
fc, approaching the MIT description. 

Remaining in the framework of non-relativistic hyper- 
radial potentials it seems reasonable to restrict the de- 
tailed analysis of all existing experimental data on single 
7-production in pp collisions to the case of g(kx) originat- 
ing from © with k = 0.1. In Fig. [5] we show the cross 
section data from 8 experiments [32] varying both in y/s 
and in the observed px region. The mean transverse mo- 
menta of the partons, necessary for a perfect description 



Pr(GeV) 

Fig. 4. The mean transverse momentum (kr) needed for the 
description of the PHENIX pp differential cross section as a 
function of pr of the outgoing photon, using the four different 
partonic transverse momentum distributions shown in Fig. [3] 




0.4 0.6 

<k T > (GeV) 



Fig. 5. The variance (fcy) as a function of {kr} for the consid- 
ered distributions g(kr)- 



of these data is plotted in Fig. [JJ We clearly see that only 
the region of small values of the ratio ^= requires rela- 
tively large (kx) for the data description. For ^ > 0.05 
the necessary mean transverse momenta lie in the interval 
[0,0.3] GeV which is in accordance with proton's struc- 
ture. 



5 Concluding remarks 

In this work we have investigated the influence of three- 
body forces in the transverse momentum distribution of 
partons inside the proton. Using a class of hyper-radial 
potentials @ we have determined the corresponding sin- 
gle parton transverse momentum distribution within a 
non-relativistic treatment having as constraints the accu- 
rate description of the proton's ground and first excited 
state energy. The charge radius of the proton turns out 
to be almost the same (~ 0.6 /m) for all potentials in 
the considered class. The obtained transverse momentum 
distributions have been incorporated in a phenomenolog- 
ical scheme, based on perturbative QCD, for the cross 
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Pt (GeV) 

Fig. 6. The differential cross section for the production of a 
single photon with transverse momentum pr in pp collisions 
as measured in various experiments. The dashed lines display 
results obtained using the resummation technique 21,22,23 . 
The results obtained through fitting of the data using the par- 
tonic transverse momentum distribution of Fig. [3]witli k = 0.1 
and variable (kr) practically coincide with the experimental 
data. 
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Fig. 7. The mean transverse momentum (fcr) needed for the 
description of all available differential cross section data for 
prompt photon production in pp collisions as a function of the 
Pt 1 1 \fs of the outgoing photon, using the k = 0.1 partonic 
transverse momentum distribution shown in Fig. [3] 



section calculation of prompt photon production in pp- 
collisions. In particular, using the associated mean trans- 
verse momentum (fcy) as a free parameter, we have fitted 
the PHENIX cross section in a wide region of the trans- 
verse momentum of the produced photon. Within this 
treatment the smoothest distribution of the (fcr)-values, 
necessary for a successful description of the data, is found 
using the transverse momentum distribution correspond- 
ing to the potential ^ with k = 0.1. This distribution has 
been also used for the description of all available experi- 
mental data for prompt photon production in pp collisions. 
The (fcT)-spectrum leading to a perfect description of all 
available experimental data is found to be restricted in 
the range [0.0,1.0] GeV. When compared with the anal- 
ysis found in the literature concerning the description of 
the same data using a Gaussian transverse momentum dis- 
tribution the results found here possess two advantages: 
(i) the interval of the necessary (fcj^-values is clearly nar- 
rower and (ii) it is displaced to smaller values which are 
closer to the geometrical characteristics of the proton ac- 
cording to Hcisenberg uncertainty relation. In a similar 
treatment in [17, using the relativistic MIT bag model we 
have obtained an even shorter interval of (fcr)-values ap- 
proaching the (fcx) ~ region. Suitably defined measures 
for the quality of the behavior of the function (hrjipr) in 
the different experiments allow for a comparison between 
the various models and lead to the following conclusions: 

— In general quark confinement leaves imprint in the 
cross-sections for prompt photon production through 
the partonic transverse momentum distribution. 

— Relativistic effects are important as dictated by the 
results found in |17] using for the confinement descrip- 
tion the MIT bag model. 

— Further study is needed in order to clarify to what ex- 
tent the exact form of asymptotic freedom (the shape 
of the inter-quark potential for small distances) is also 
influencing the quality of the description of experimen- 
tal data within our approach. 

— Finally it turns out that three-body forces, included in 
the present approach but not in the MIT bag model, 
are also important for an efficient description of par- 
tonic transverse momentum effects inside the proton. 

Thus it is interesting to extend the present work by inves- 
tigating the partonic transverse momentum distribution in 
a MIT bag model with interacting partons where three- 
body forces are also included. According to the findings 
of the present work such a model should lead to a further 
improvement of the description of the proton transverse 
momentum structure. 
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